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ABSTRACT 

High-quality averaged radio profiles of some pulsars exhibit double, highly symmetric 
features both in emission and absorption. It is shown that both types of features are 
produced by a split-fan beam of extraordinary-mode curvature radiation (CR) that 
is emitted/absorbed by radially-extended streams of magnctospheric plasma. With 
no emissivity in the plane of the stream, such a beam produces bifurcated emission 
components (BFCs) when our line of sight passes through the plane. A distinct ex- 
ample of double component created in that way is present in averaged profile of a 5 
ms pulsar J1012+5307. We show that the component can indeed be very well fitted 
by the textbook formula for the non-coherent beam of curvature radiation in the po- 
larisation state that is orthogonal to the plane of electron trajectory. The observed 
width of the BFC decreases with increasing frequency at the rate that confirms the 
curvature origin. Likewise, the double absorption features (double notches) are pro- 
duced by the same beam of the extraordinary-mode CR, when it is eclipsed by thin 
plasma streams. The intrinsic property of CR to create bifurcated fan beams explains 
the double features in terms of very natural geometry and implies the curvature origin 
of pulsar radio emission. Similarly, "double conal" profiles of class D are due to a cut 
through a wider stream with finite extent in magnetic azimuth. Therefore, their width 
reacts very slowly to changes of viewing geometry due to the geodetic precession. 
The stream-cut interpretation implies highly nonorthodox origin of both the famous 
S-swing of polarisation angle, and the low-frequency pulse broadening in D profiles. 
Azimuthal structure of polarisation modes in the CR beam allows us to understand 
the polarised 'multiple imaging' and the edge depolarisation of pulsar profiles. 



Key words: pulsars: general - pulsars: individual: JlOl^ 
B0525+21 B1913+16 - Radiation mechanisms: non-thermal. 
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1 INTRODUCTION 

Double 'absorption' features in radio-pulse profiles were first 
identified in radio data independently by Rankin & Rath- 
nasree (1997; B1929+10) and Navarro, Manchester, Sandhu, 
et al. (1997; J0437-4715). McLaughlin & Rankin (2004) dis- 
covered the double notches in the leading wing of the main 
pulse of B0950+08. 

Navarro et al. have noticed that the feature must be 
'intrinsic to the emission mechanism' because it becomes 
narrower at larger observation frequency v b B . However, the 
initial interpretive efforts of theorists did not follow that sug- 
gestion. Wright (2004) interpreted the features in terms of 
altitude-dependent special-relativistic effects, but was forced 
to assume unlikely emission geometry and postulated an 
opaque absorber of unknown origin that corotates at/near 
the light-cylinder. This work, however, pioneered the im- 
portant idea that a single entity must be responsible for 
both notches. Dyks, Frackowiak, Slowikowska, et al. (2005) 



considered the neutron star itself (embedded in an opaque 
plasma cloud) as the absorb er/eclipser and their emission 
region neatly followed the geometry of magnetic field lines. 
However, to make this far more natural geometry really 
work, the radio emission had to be directed inward, toward 
the neutron star. The model has achieved some agreement 
with the data, but the symmetry of double notches, as well 
as their frequency evolution remained unsolved. 

Considerable progress in the ability to interpret the 
data was done by Dyks, Rudak & Rankin (2007; hereafter 
DRR), who interpreted the notches as a direct imprint of 
the hollow beam intrinsic to the radiation emitted by elec- 
trons accelerated parallel to their velocity. This work was 
the first to interpret the double features in terms of micro- 
physical beam intrinsic to a specific radiation mechanism. 
The model managed to ensure symmetric and double shape 
of the notches whenever they were not washed out by spa- 
tial extent of the emitter/absorber. The separation A n of 
the notches (as observed at that time) was consistent with 
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Figure 1. Pulse profiles of PSR J 0437— 4715 at three frequencies. 
The double features at rf> ~ 51° and 70° make it possible to 
align the profiles in absolute way. Note the lack of "radius-to- 
frequency" mapping. Zero levels and flux units are arbitrary and 
different for each profile. Data courtesy: R. Manchester. 



the inverse-Compton-like version of the parallel acceleration 
maser: ^ h s ~ 7 2 ^pi A^ 2 '. Since the mechanism is intrin- 
sically broad-band, one could naturally understood the lack 
of radius-to- frequency mapping in J0437— 4715 (Fig. [TJ. 

DRR noticed that double notches of J0437— 4715 are 
located in a trailing wing of a bifurcated emission compo- 
nent (BFC), visible at (j> ~ 51° in Fig. [1] This component 
was interpreted in terms of the same elementary emission 
beam that produces the double notches. Throughout this 
paper we maintain this view that the apparently unrelated 
bifurcated emission components and double notches have 
common origin. The notches are regarded as a negative im- 
age of the elementary emission beam of the same radiative 
process that produces the bifurcated components. 

In spite of the overall success of DRR's model, the par- 
allel acceleration beam was unable to explain the observed 
large depth of double notches. To produce the notches, a 
small part of the emitter had to be radio silent or a lo- 
calised absorber of unknown origin was needed. Neither it 
was possible to decipher the real macroscopic geometry of 
the emitter/absorber system. 

In this paper we solve these problems and we identify 
both the topology of the elementary emission beam as well 
as the general geometry of the eclipsing phenomenon. This 
is done by making physical fits to a bifurcated emission com- 
ponent in the averaged profile of PSR J1012+5307 (Section 
I3.2|l . as well as by 3D simulations of double notches (Section 
2| . Uncovering of the nature of double features is equivalent 
to the identification of the long-sought radio emission mech- 
anism of pulsars, and results in important consequences for 
several problems of radio pulsar astronomy (Section |SJ). 
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Figure 2. Pulse profiles of J1012+5307 at two frequencies. They 
are phase-aligned to match the pronounced bifurcated component 
(BFC) near <f> ~ 71°. This also aligns the central minimum and 
the trailingmost component of the main pulse (MP). Note the 
other weak BFC located ~50° ahead of the MP in the 0.82 GHz 
profile (around </> ~ 245°). 

2 THE BIFURCATED COMPONENT OF 
J1012+5307 

The 5.25-ms pulsar J1012+5307 was discovered by Nicas- 
tro et al. (1995). Its polarisation properties were studied by 
Xilouris et al. (1998) and Stairs et al. (1999). Kramer et 
al. (1999) investigated the behaviour of this pulsar across 
the radio frequency spectrum, and it was the high-quality 
1.4-GHz profile presented in that paper that has drawn our 
attention to this object. 

Fig- EH presents the averaged pulse profiles of 
J1012+5307 observed with the Green Bank Telescope at 
0.82 (bottom) and 1.4 GHz (top). The bandwidth was 64 
MHz at both frequencies. The total integration time of ~15 
hours (~10 7 single pulses) was accumulated in the period 
between July 2004 and March 2007. A pronounced, highly 
symmetric BFC can be seen near the phase <j> = 70° (~50° 
ahead of the interpulse). Another bifurcated component (a 
weak one) is present some 50° ahead of the main pulse. 

The outstanding signal-to-noise ratio of the data, and 
the strength of the bright double component make it possible 
to test quality of various empirical or physical fits. In Fig. [3] 
we present the traditional decomposition of the BFC into 
a sum of two gaussians. The phase interval of the fit was 
limited to the central part between the vertical dashed lines. 
One can see that modelling of double features with gaussians 
is not a good idea. The problem is that any sum of two 
symmetrical functions tends to have the inner wings less 
steep than the outer wings, which is just opposite to what 
is observed Q 

1 The gaussian fitting is also problematic for single asymmet- 
ric components, as was rightfully emphasized by Weltevrede & 
Johnston (2008a). Still, the authors decided to use symmetric 
functions. 
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Figure 3. A best fit of a sum of gaussians to the BFC of 
J1012+5307 at 0.82 GHz. The vertical dashed lines (also in the 
next figures) mark the phase interval in which the fit was per- 
formed. In this interval x 2 /dof ~ 14. 




i i I [ i 

-20 -10 10 20 

pulse phase [°] 



Figure 4. Best fit solution for the off-axis cut through the 
frequency-integrated parallel-acceleration beam (eq. [T}- Within 
the marked phase interval x 2 /dof = 2.2. 



3 PHYSICAL FITS TO A RADIO PULSE 
COMPONENT 

Because of astonishing degeneracy of Nature the BFC com- 
ponent can be fairly well fitted with two different physi- 
cal models. To expose the deceiving character of the prob- 
lem, we intentionally start the analysis with the incorrect, 
parallel-acceleration model (hereafter PAC model). 



3.1 The wrong idea: parallel acceleration beam 

At 0.82 GHz, the central minimum in the BFC component 
of J1012+5307 reaches a large depth of nearly 70%, which in 
the PAC scenario could only be interpreted as emission from 
a very localised region. It is therefore assumed that the BFC 
is free from any spatial convolution effects and we model it 
directly with the hollow parallel-acceleration beam. 

We use the textbook dipolar radiation pattern for non- 
coherent radiation emitted by relativistic charges acceler- 
ated parallel to their velocity (eg. Rybicki & Lightman 
1979). The shape of the beam (frequency-integrated power 
per unit solid angle) is 



sin 2 6 



Aire 3 



(l-pcosoy 



(i) 



where is the angle between the beam axis (electron veloc- 
ity) and the line of sight, /3 = (1 — l/'y 2 ) 1 ^ 2 is the electron 
velocity in units of the speed of light c, whereas q and ay 
are the charge and acceleration of the electron. 

This beam is inclined at some angle ctb with respect to 
the rotation axis and our line of sight passes through it at 
some 'beam impact angle' fit = C — a b- The beam inten- 
sity ([l]) is related to the observed pulse longitude (phase) 
through the cosine rule for sides of the spherical triangle 
(a,v,H): 

(2) 



9 — cos(</> — cpo) sin otb sin ( + cos Q(, cos 



cos 

where <j>o is the phase of the central minimum in the bifur- 
cated component, Q, is the angular rotation velocity of the 



pulsar, v is the velocity of the emitting electron(s), and n is 
the unit vector of our line of sight. 

The zero-flux level is blindly fixed at the lowest place 
in the profile. Since the inclination of the beam axis merely 
rescales the feature by a factor 1 / sin at we assume ab = 90° . 
Thus, the fitted function F((f>) = CP\\ (cos 6(<f))) has 4 param- 
eters: the normalisation constant C, the phase location 0o of 
the double component, the viewing angle £ (or the beam im- 
pact angle f3b), and the electron Lorentz factor 7. The fitting 
was done with the Levenberg-Marquardt method. To make 
sure that the best fit corresponds to the global minimum, 
the fitting was repeated for many sets of initial parameters. 
These were sampled quasi-randomly from wide intervals, us- 
ing the maximum-avoidance algorithm of Sobol', Antonov, 
and Saleev (Press, et al. 1992). 

The best-fit solution for the 0.82 GHz profile is shown 
in Fig. [4] The fit was performed for the phase interval which 
is marked with the vertical dashed lines and has a width of 
A<f) — 1.8Abf c , where Abf c is the separation of maxima of 
BFC. Within this interval xVdtf stays close to unity, and 
it is equal to 2.2 for the specific case shown in Fig. [4] The 
beam impact angle is f3b = 1.546 ± 0.010°, and the elec- 
tron Lorentz factor 7 = 6.364 ± 0.012. For a non-orthogonal 
ctb 7^ 90° this translates to a somewhat larger value of 
7 ~ 6.4/ sin a;,. The maxima are separated by the phase 
interval Abfc = 8.54 ± 0.020° and the opening angle of the 
hollow cone is 2# max = 9.08±0.030° (or 9°/ sin ab in the case 
of ab 7^ 0). The errors include only the statistical la uncer- 
tainty and are underestimated due to the unknown zero-flux 
level. The zero-level is likely to be wrong by few percent of 
the flux observed at the maxima of the BFC component 
(/max)- This is suggested by the flux difference of ~0.01/ma X 
observed at two lowest locations in the 1.4-GHz profile (at 
(f> ~ 30° and 200° in Fig. [5J. The corresponding system- 
atic errors were estimated by repeating the fit for various 
levels of the zero flux. A shift by 5% percent of /Jfax typ- 
ically changed the values of 7, # max and Abf c by 1 — 2%. 
The beam impact angle f3b (and the normalisation constant 
C) were more sensitive and varied by roughly the same fac- 
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Figure 5. Best fit (x 2 /dof ~ 12, thin line) of the parallel ac- 
celeration beam to the BFC at 1.4 GHz. The asymmetric thick 
line presents a composition of two separate fits to the left- and 
right-hand part of the BFC. In this composed case x 2 /dof ~ 1. 



tor (5%). The actual errors of the parameters are therefore 
dominated by the systematic effects of unknown zero level 
and have the magnitude of few percent. Below, the fitted 
parameters will be given without errors. 

In the central parts of the BFC, where bulk of the emis- 
sion is received, the parallel acceleration beam reproduces 
the data with very high accuracy. This is astonishing, be- 
cause the fitted function ignores all spatial and spectral con- 
volution effects (note that frequency-integrated emissivity 
is fitted to a component observed with fairly narrow band- 
width of 64/820 MHz = 8%). The fitted curve also ignores 
all deformations that could possibly arise in the amplifica- 
tion process. 

This apparent coincidence is very deceiving and can eas- 
ily mislead one to believe in the parallel-acceleration sce- 
nario. However, it merely tells us that the component can 
be well modelled with the function of type 6 2 /(l — j3 cos 6) 3 , 
with the exponent g weakly constrained, because it mostly 
determines the shape of the outer wings. The PAC is not 
the only mechanism described by such function. So the re- 
sult does not prove the parallel acceleration origin. It rather 
suggests that the amplification process preserves the shape 
of the non-coherent radiation beam. 

The fit is worse at 1.4 GHz (thin line in Fig. [5]), 
where the BFC is more asymmetric. However, if the left- 
and right-hand half of the BFC is fitted separately (thick 
line), the fit quality is similar (x 2 /dof ~ 1 again within 
Acb ~ 1.8Abf c ). The parameters obtained with the single- 
step fit are: j3 b = 1.7°, 7 = 7.48, 2<9 max = 7.7°, A bfc = 6.9°. 
The result again looks misleadingly encouraging: the beam 
impact angle is, within the few percent errors, the same as 
at 0.82 GHz. Thus, using the parameters obtained from the 
fit at 0.82-GHz, it is enough to change only a single param- 
eter (the Lorentz factor 7) to reproduce both the relative 
depth of the central minimum and the separation of max- 
ima at 1.4 GHz. However, the ratio of 71.4/70.82 = 1.17, 
that is required to achieve this, is smaller than expected for 
the inverse-Compton-like process: the inverse-Compton re- 
lation !/ b s oc 7 2 implies 71.4/70.82 = 1.3. In other words, 
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Figure 6. Merging-rate index a in the relations: A^fc oc v £ 
and A n oc k'T", where Ai-,f c is the separation of maxima in bifur- 
cated components (BFC), whereas A n is the separation of min- 
ima of double notches (DN). The vertical lines mark the values 
of a for the curvature radiation (dotted) and inverse-Compton- 
like emission due to the parallel acceleration (dashed). The errors 
are statistical ltr only. Note that the only well-resolved and well- 
modelled case is the BFC of J1012+5307. 



the inverse-Compton scenario requires that the beam size 
should decrease with increasing frequency as # ma x oc v~^ 5 
whereas the inferred rate is 6 maK oc t^g' 31 . Thus, with the 
inverse-Compton-implied ratio of 71.4/70.82 = 1.3, one can- 
not fit the data simultaneously at two frequencies. Instead, 
the separation of maxima in the BFC follows the relation 



Abfc oc (see Fig. [6} which is fully consistent with 

the low-frequency curvature emission for which the opening 
angle is i/Wx oc v~£/ 3 (eg. Jackson 1975) . 

A much more serious problem with the PAC model is 
that the intrinsic topology of the beam (axially-symmetric 
hollow cone) makes it extremely difficult to obtain double 
notches with the observed depth of >20% (Perry & Lyne 
1985; see also fig. 2c in DRR). 



3.2 The curvature radiation 

3.2.1 The classical radiation beam 

In the case of curvature radiation, the energy emitted per 
unit frequency and unit solid angle is given by: 



f]crv 



m\ + n ± 



2 2 

q ui 
3tt 2 c 



( E Y k 2 Kl(y)+Z Kl(y) sin 2 ^ 

\Cj 13 3 



where 

C = l/ 7 2 +^ 2 , 

^P c 3/2 

V= 3c^ 



and 



(3) 
(4) 

(5) 
(6) 



(eg. Jackson 1975). The symbol p denotes the radius of cur- 
vature of electron trajectory in our reference frame, ui = 
2m/ b s , c is the speed of light, and K's are modified Bessel 
functions. The angle ip is measured between our line of sight 
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Figure 7. a) Transverse crossection through the fan beam of 
low- frequency curvature radiation (eq. |4}. The angle ip ls mea- 
sured from the plane ip = of electron trajectory (the beam does 
not have axial symmetry). The total emissivity r\ (thin solid) con- 
sists of the 1 1 -mode radiation that peaks in the trajectory plane 
(thin dot-dashed line) and the two-planar _L-mode which has the 
double-peaked crossection shown with the thick solid line. The 
thick dot-dashed lines (marked ||A, ||B, and ||C) show possible 
profiles of partially-attenuated ||-mode radiation, b) The polari- 
sation degree L/I and the ratio fjx/^y f° r the case with no at- 
tenuation of the ||-mode. Note the intrinsic decrease of L/I in the 
outer wings of the fan beam (see discussion in Sections l5land ltT2)l . 



and the plane of electron trajectory^ The "well-known" 
shape of the beam is shown in Fig. [7] The beam can be 
decomposed into two, orthogonally polarised sub-parts: the 
one which is polarised parallel to the projection of electron 
trajectory on the sky (r/\\) and the other one (r/±), with the 
polarisation orthogonal to the plane of the electron trajec- 
tory (hereafter ET plane) 



2 The angle ip should not be mistaken with the polar angle 8 
measured from the electron velocity or B. 

3 There can be a big difference between the shape of electron 
trajectory in the observer's frame and the shape of the i?-ficld line 
along which the electron propagates. However, it is not important 
for the considerations of the present paper, and will be ignored 
in what follows. Thus, the ET plane can be considered to be the 



As shown in Fig. [7] all the radiation in the plane of 
electron trajectory (ip — 0) is polarised parallel to the ob- 
served projection of the B-field line and is therefore likely to 
have problems with leaving the magnetospheric plasma. As 
can be seen in eq. @, the remaining part of the beam (the 
orthogonal mode, hereafter _L-mode) has the central mini- 
mum of exactly the same shape (as a function of ip) as the 
PAC beam has as a function of 9 (near the B-field line we 
have Py oc sin 2 8 for the PAC model, and tj± oc sin 2 ip for 
the _L-mode curvature radiation). The double shape of the 
orthogonal-mode beam makes it a suitable tool to model the 
bifurcated components. 



3.2.2 Fitting the curvature beam 

We assume that the emitter has a form of a thin and elon- 
gated plasma stream that emits the curvature radiation 
mainly in the extraordinary (orthogonal) mode. At each al- 
titude within some finite but non-negligible range of Ar, a 
broad range of frequencies is emitted. Therefore, the beam 
observed at a fixed frequency has a fan-like shape and sub- 
tends a range of magnetic colatitudes 9 m . Since there is no 
_L-mode emission in the plane of the stream, the BFC com- 
ponent is naturally produced when our line of sight passes 
through the plane of the stream. 

We proceed with the fitting as in the parallel accelera- 
tion case, ie. we use the frequency-integrated power emitted 
per steradian. Integrating eq. Q, and dividing by 2irp/(cp3) 
one obtains: 



IcTV, II -^crv,^ 

7q 2 c ( 2/3 7 x 1/2 



2567T/9 2 \ cosi/j 

/, r, ,%-s/2 / ■ , ,s 5 (3 cos ip sin 2 tp 

x l_/3cosy< J/2 UWO + Tjn — I — r 
V 14 1-/3 cos tp 



(7) 



,(8) 



(eg. Konopinski 1981). The function h(ip) describes the un- 
known contribution of the ordinary (parallel) mode and is 
equal to 1 in the absence of absorption. For a pure orthog- 
onal mode, ie. for I CIV — J crv .i, we have h(ip) — 0. 

To record the bifurcated component, the line of sight 
must pass through the plane of electron trajectory. In the 
case of the pure orthogonal mode this implies zero flux at 
the center of BFC. In reality, however, there are several fac- 
tors that are likely to raise the central minimum. These in- 
clude a small (but non-zero) thickness of the plasma stream, 
partial admixture of the ||-mode, and non-planarity of the 
trajectory. To be able to fit the central part of the BFC 
we artificially raise the level of the central minimum in two 
ways: 1) h(ip) is considered as a ip- independent normalisa- 
tion parameter for the parallel mode (h(ip) = ho < 1); 2) we 
redefine the meaning of the angle ip. 

For numerical simplicity, we assume that the plane of 
the stream is meridional, and the general geometry of the 
problem is quasi-orthogonal: ctb = 90°, C, ~ 90°. In such a 
case ip simply becomes 



ip ~ 4> - cj> , 



(9) 



plane of the S-field line along which the electron propagates in 
the corotating frame. 
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Figure 8. Best fit of the frequency-integrated and _L-mode- 
dominated curvature beam (eq. HOC to the 0.82-GHz BFC of 
J1012+5307. Note the much wider interval of the fit, in which 
X 2 /dof ~ 3.6. 



where <j> is the observed pulse phase and 4>o is the phase at 
which the center of the BFC is observed. So the situation is 
different than in the case of eq. ([2]). However, it is convenient 
to mimic the non zero flux at the center of the profile by 
assuming cos tp ~ cos 8 as given by eq. ((2} . This equation 
reduces to eq. ([9]) whenever at — C <S 4> — 4>o, ie. everywhere 
except near the central minimum, since in our case Ob — f — 
1.65°. The value of ( fitted in that way has no physical 
or geometrical meaning. It is an arbitrary parameter that 
allows to mimic the unknown contribution of the ordinary- 
mode radiation (or other effects). 

Thus, the fitted function was either given by eq. ([8]) 
with ip = cf> — c/>o, or by the following formula: 



A 



3 1 ' 2 ip sin 2 ip 



(1-13 cos </>) 7 / 2 ^ 

with cosi/' = cos# given by eq. ([2]). Both methods gave 
similar parameter values and had comparable x 2 /dof, with 
eq. (fT0)l doing only slightly better than ([8| . 

Fig. © shows the fit of eq. (JTOJ to the 0.82 MHz profile. 
The agreement is only slightly worse than in the case of the 
parallel acceleration beam (x 2 /dof = 3.6) but holds within 
a much larger interval of phase (marked with the vertical 
dashed lines). A fit at 1.4 GHz (Fig. [9]) is worse due to the 
BFC's asymmetry and likely due to increased contribution 
of the ||-mode. 

The extraordinary curvature beam works better in the 
outer wings because of the much smaller exponent (g = 7/2) 
in the denominator of eq. (imp (cf. eq. [T]in which g = 5). 
[Note: Because in the outer wings we have 1 — /3 cos tj) ~ 
i/> 2 /2, the wings have the same shape of (1 — f3 cos tp)~ 5 ^ 2 
for both the parallel and orthogonal polarisation state (see 
eq. [8} • Therefore, we propose that single pulse components 
be fitted with the function of type: 

!(</>) = J 4(l-/3cos(0-0 o ))" 5/2 (11) 

-5/2 



" A I^ + W-«J . (12) 
where /3 = (1 — I/7 2 ) 1 / 2 , and 7, A, and 0o are the parameters 
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Figure 9. The curvature-radiation fit at 1.4 GHz. 



of the fit. The value of 7, which determines the steepness of 
wings should be different for the leading and trailing half of 
an asymmetric component. Such function is very successful 
in reproducing the wings of profiles and allows to avoid the 
false components that usually come out in the procedure of 
gaussian fitting (see fig. 2 in Weltevrede & Johnston 2008a 
and the comments therein).] 

However, the agreement visible in Fig.[5]is not ideal and 
the fitted function again ignores the issue of limited band- 
width as well as all possible convolution effects. Therefore, 
the quality of the fits shown in Figs. [3] and [S] if considered 
alone, does not allow us to firmly discriminate between the 
PAC maser and the curvature radiation. What indicates the 
curvature origin of the beam, is that the fitted separation of 
peaks (8.30° and 6.87° at 0.82 and 1.4 GHz) is consistent 
with the relation Abf c oc h> ohs ' . Such relation is expected 
if the observed frequency is smaller than, or equal to the 
characteristic frequency u crv of the curvature spectrum: 



7 



= 0.7 MHz 



(7/10) 3 
p/(10 7 cn 



(13) 



^crv = | £ ^=7-10 9 H Z - , 

4-7T p p/(lcm) 

In both these cases (v b s <S ^crv, or i/ Q b s — v C iv) the open- 
ing angle of the beam is i/> max ~ [3c/(27rp^ b s )] 1 ^ 3 , which 
follows directly from the properties of the fixed- frequency 
radiation pattern (eq. [4] eg. Jackson 1975). The exponent of 
1/3 matches the behaviour of the BFC of J1012+5307: the 
observed Abf c can be well-fitted simultaneously at two fre- 
quencies, with the ratio of 71.4/70.82 fixed at the predicted 
value of (1400/820) 1/3 = 1.195. However, the exponent of 
1/3 is not ubiquitous among pulsars (see Fig. |SJ). An addi- 
tional criterion is evidently needed to discriminate the mod- 
els. 



4 DEPTH OF NOTCHES AS THE DECISIVE 
DISCRIMINATOR 

As shown in DRR the PAC maser is unable to naturally 
produce double notches with the observed depth of > 20% 
(Perry & Lyne 1985). This is because the PAC beam has 
the hollow cone shape, and is axially symmetric around the 
local electron velocity (around the local B). Therefore, a 
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Figure 11. The beam of the previous figure is carried by the elec- 
tron along a magnetic field line, which produces the two-planar, 
split-fan beam. The fixed-intensity contours of the instantaneous 
beam are shown at three locations on a magnetic field line. 



Figure 10. Instantaneous beam of _L-mode curvature radiation 
(power per unit frequency per unit solid angle) as a function of the 
angle ip, measured from the plane of a B-ficld line, and the angle 
fi measured in the orthogonal way. The vector of electron velocity 
(not shown) is anchored at the origin of the reference frame and 
points vertically upward, along the local magnetic field. There is 
no radiation in the plane of the electron trajectory, which is the 
plane of -0 = 0. 

ring-shaped part of the emitter contributes to the flux that 
is observed at the minima of notches and makes them very 
shallow. For a two-dimensional emitter that extends later- 
ally at a fixed altitude, this 'light-polluting' ring has the 
crossection which is marked with letters 'b' in fig. 3 of DRR. 

In the case of the _L-mode curvature radiation, the beam 
has the mirror symmetry with respect to the plane of the 
local B-field line, and there is no (or little) emission within 
the entire plane (see Figs. llOl and lll[) . Each point of emission 
region is therefore mostly emitting in only two directions, 
which make small angles i/> ma x with respect to the plane of 
the £?-field line (note that in the case of the low-frequency 
curvature radiation ip msx S> I/7). Therefore, instead of the 
'polluting ring', there are only two spots in the emission 
region (on both sides of the absorber) that contribute to the 
flux observed at the minima of double notches. The minima 
are therefore much less contaminated by the nearby parts of 
the emitter than in the PAC case. Moreover, the eclipsing 
object does not longer need to be compact. Overlying opaque 
streams (Fig. I12[) . or elongated disruptions in the emitter, 
can produce deep notches without affecting the flux at the 
center of the notches. In the curvature scenario deep double 
notches are therefore produced easily and naturallyQ 

Fig. [11] schematically shows how the split-fan beam 
due to a thin stream is created: when electrons move along 
the curved magnetic field lines they carry-along the two- 

4 After many failed numerical attempts to reproduce the ob- 
served depth of notches (using the axially-symmetric beam of 
PAC), DRR concluded: "The large depth of double notches could 
easily be produced if the coherent radio emission had occured at 
two small angles (I/7) with respect to the plane of B-field lines." 




Figure 12. Natural geometry that produces deep double notches. 
An opaque stream (black arc) is located above a two-dimensional 
emission region. The solid lines present magnetic field lines. Two 
dashed lines trace the eclipsed parts of the emitter. 

directional pattern of Fig. [TO] The resulting beam has 
therefore completely different structure (split fan) than the 
parallel-acceleration beam (hollow cone). It is precisely this 
topological property that really makes the crucial difference, 
and allows us to recognize which mechanism is the correct 
one. 

The ease with which the _L-mode curvature beam pro- 
duces deep double notches, was verified with a 3D numerical 
simulation. The code assumes that the emitter has the shape 
of a surface formed by magnetic field lines with fixed foot- 
print parameter s. The footprint parameter is arbitrarily 
fixed at s = 1 (last open field lines). To produce the om- 
nipresent 'pedestal' emission, the curvature radiation con- 
tinues up to a quite high altitude of (2/3)i?i c Q 

The most natural way to generate the notches is shown 

5 The pedestal can also be created by low-altitude emission from 
B-field lines with s > 1 (closed field line region) or even from 
fixed-altitude region involving a range of s. 
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in Fig. [T2]and involves a thin opaque stream of plasma above 
our emission region. We call it most natural because we know 
that such streams exist and produce the BFC. Note that be- 
cause the notches are observed in averaged profiles, the ob- 
scuring stream has to be fairly permanent to be present in 
the same pulses as the emission which is absorbed. However, 
to simplify the calculations we preferred to model an elon- 
gated fissure in the emitter, instead of the obscuring stream. 
Thus, we assumed that a narrow part of the emitter (a wedge 
with a fixed angular width of A0 m = 4.1° in the magnetic 
azimuth) is not emitting. This is a safe assumption, since the 
ability to produce deep notches does not depend on whether 
they are produced by an absorber above, or the hole within 
the emission region (DRR). 

The simulations have confirmed that deep double 
notches appear naturally for the _L-mode curvature beam. 
A sample result for three different viewing angles is shown 
in Fig. 1131 It also becomes clear now why the notches can 
only be observed in very highly polarised pedestal emission: 
when the ||-mode is present, it fills-in the central minimum 
of the double _L-mode beam and the notches cannot appear. 

Based on Sections 3 and 4 we learn that in the case 
of curvature radiation, both the double absorption features 
(double notches) of several pulsars as well as the bifurcated 
emission components of J1012+5307 and J0437— 4715 can be 
understood in terms of natural geometry that involves the 
absorption/emission by plasma streams in pulsar magnetic 
field. We conclude that it is the curvature radiation, not the 
radiation due to the parallel-acceleration, that is observed 
as the coherent radio emission from pulsars. 




-40 



-20 



20 



40 



Figure 13. Deep double notches simulated with a 3D numerical 
code. The background (pedestal) emission was created by the 
surface of last open field lines. A wedge-shaped part of the surface 
was not emitting to produce the notches. The wedge was centered 
at the magnetic meridian 4> m = —20°, and had the azimuthal 
width of Aifim = 4.11°. The curves correspond to £ = 64.3°, 
51.5°, and 38.6° (top to bottom). The magnetic dipolar field was 
tilted at the angle a = 90° . 



5 POLARISATION STRUCTURE OF THE 
CURVATURE BEAM 

To consider some implications of our findings, we first dis- 
cuss the polarisation structure of the two-modal curvature 
beam. 

It is convenient to imagine the 'elementary' beam as 
the fan beam due to the passage of an electron along a mag- 
netic field line (or due to a thin stream flowing along a nar- 
row bunch of i?-field lines) . Based on the fits of Section [3] 
it is reasonable to assume that the _L-mode beam is am- 
plified isotropically and preserves the original (non-coherent 
shape). The ||-mode beam is likely to be attenuated in the 
central parts (near tp ~ 0), where the radiation propagates 
through the largest plasma density. The outer parts of the 
||-beam are emitted from the outer boundary of the stream 
into the low-density region surrounding the stream. Guided 
by the observed edge depolarisation of profiles (see Section 
I6.2|l we assume that in the outer wings of the beam, the ||- 
mode radiation is amplified with the same efficiency as the 
_L-mode and is not attenuated. 

The polarisation-structure of such a toy beam is shown 
in Fig. [7^, with thick lines (solid for _L-mode, dot-dashed 
for ||-mode). Three cases of ||-mode attenuation are shown, 
marked with ||A, ||B, and ||C. Since the ±-mode has zero 
emissivity in the ET plane, even a small contribution of 
||-mode (shown as the case |]B) can dominate the central 
parts of the total beam. It is also possible to have the ||- 
mode fully absorbed in the beam center (case ||C). Such 
half-phenomenological curvature-radiation beam has several 



intrinsic properties that make it notably successful in inter- 
preting enigmatic pulsar phenomena^ 



6 APPLICATION TO SELECTED PULSAR 
PROBLEMS 

6.1 Polarised multiple imaging 

The inferred polar carousel of many pulsars exhibits inter- 
esting structure in magnetic azimuth (f> m : subbeams of op- 
posite orthogonal polarisation are distributed alternatingly 
along the polar ring of emission. This pattern produces sev- 
eral enigmatic phenomena, such as 1) the "modal parity at 
the outside edges of polar beam" (Rankin & Ramachandran 
2003), 2) two orthogonally-polarised, out-of-phase drift pat- 
terns, 3) jumps in drift phase (Edwards 2004) . The latter two 
effects were described by Edwards, Stappers & van Leeuwen 
(2003) as 'multiple imaging' possibly caused by refraction, 
birefringence or aberration/retardation effects. 

Without intending to model the phenomena in detail, 
we want to show that the low-frequency curvature radia- 
tion beam has several intrinsic properties that are consis- 
tent with the observations. First, the multiple imaging is 

6 Progress in understanding of the polarisation structure could 
be achieved by decomposing the BFC of J1012+5307 into orthog- 
onal polarisation modes. Unfortunately, the polarisation data for 
J1012+5307 are not available to us at the present moment. 
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automatically implied by the non-zero angular width of the 
low-frequency beam of curvature radiation: 



0.62 



/ 3c 



V 2-wpv ohs 



1/3 



0.40°(p7^)" 1/3 , 



(14) 



where p7 = p/(10 7 cm), and V9 — f bs/(l GHz). One can 
see that for the typical radius of curvature of 10 7 cm, the 
intrinsic width effects become important/noticeable right in 
the lower part of the standard frequency band of radio as- 
tronomy (fobs < 1 GHz). Any intrinsic broadening of com- 
ponents that is smaller than ~1° would be hard to notice in 
pulse profiles. Thus, by transiting from sub-GHz to super- 
GHz radio band one goes from profiles dominated by in- 
trinsic shape of the curvature-radiation beam (with the ra- 
diation emitted at noticeable angles with respect to B) to 
profiles dominated by the spatial distribution of emission 
region (with the emission tightly bound to the direction of 
B). At fixed fobs the intrinsic broadening depends only on 
the curvature radius p and the dependence is rather weak 
(cc p 1 ^ 3 ). This explains why the scale of double features is 
so similar in so different objects like the millisecond pulsars 
(MSPs) and normal pulsars. On the other hand, p is a bit 
smaller in MSPs and the strongly flaring B-field lines make 
it possible to observe the double emission components well- 
isolated from the main pulse. Therefore, the double emis- 
sion features are most easily noticed in MSPs. The scale of 
double absorption features (notches) additionally depends 
on the eclipsing geometry and will be analysed elsewhere. 
We conclude that the multiple imaging effect is naturally 
guaranteed by the non-zero width of the curvature beam at 
low fobs- No refraction, birefringence etc. is needed in the 
zeroth-order approach to understand the multiple imaging. 

6.1.1 A method to determine the curvature of electron 
trajectory in the emission region 

At a fixed and low observation frequency (f bs ^ ^civ), the 
opening angle VWx of curvature emission depends only on 
the curvature radius of electron trajectory in our (observer's) 
reference frame. This opens a way to calculate the curvature 
radius p in the emission region from eq. (|14fl . However, this 
can only be achieved when full information on the global 
geometry of the system is available. Since the viewing angle £ 
can be different than 90° , the observed separation of maxima 
in the BFC can be rescaled (increased) by the factor 1 / sin £ 
due to the well-known 'not-a-great-circle' effect. Moreover, 
if the plane of the electron trajectory is inclined at an angle 
5mda with respect to the local meridional plan^ll the value of 
Abfc is increased by additional factor 1 / cos <5 m d n • One can 
also use the factor 1/ sin 5 cu t, where 8 cut = 90° — <5 m d n is the 
angle at which we cut the ET plane0 The observed Abfc 
is therefore related to ip max through the global geometry 
parameters: 



2^n 



0.8° 



sin ( cos <5 mdn sin £ sin 5 cu t {pivg) 1 ^ 



(15) 



7 The rotational meridian at which our line of sight crossess the 
plane of electron trajectory is meant here. 

8 Thus, S cu t is the angle between the ET plane and the plane 
that is locally tangent to the surface of the cone encircled by our 
line of sight. 



Therefore: 

1.72° 



Pe 



Abfc sinCsin5cut 



(16) 



where p 6 = p/(10 6 cm). For J1012+5307 we have A b fc = 
8.61° and 7.1° at 0.82 and 1.4 GHz respectively, which im- 
plies p sin 3 £ sin 3 5 cut — 0.011 • 10 6 cm at both frequencies. 
To make p equal to 10 6 cm, either £ or <5 cu t of ~13° is re- 
quired. For moderate values of £ = <5 cu t — 45°, the value of 
p is increased by a factor of 8 up to ~0.1 • 10 6 cm. 

Note that the small value of p sin 3 £ sin 3 <5 cu t does not 
have to necessarily imply sub-stellar curvature radius. The 
main pulse of J1012+5307 has the width that corresponds 
exactly to the opening angle of low-altitude polar tube, and 
the interpulse is present. This suggests orthogonal viewing 
geometry (f ~ a ~ 90°). However, this kind of geometry si- 
multaneously necessitates small cut angles (<5 cu t ~ 0°) for 
components observed well ahead of main profile features 
(main pulse and interpulse). 

In the case of J0437-4715, we have A bfc = 3.76° , 2.89° , 
and 1.7° at 0.438, 0.66 and 1.512 GHz, respectively. This 
gives p sin 3 sin 3 5cut = 0.22, 0.3, and 0.7 • 10 6 cm, respec- 
tively. The results for J0437— 4715 are closer to what we 
expect for dipolar magnetic field lines. The actual p is prob- 
ably larger, because £ < 45° is likely. On the other hand, 
the BFC of J0437— 4715 has a very shallow central minimum 
due to some convolution effects (eg. due to the E x B drift 
of the stream). This can make Abfc appear slightly smaller. 

It is worth to digress here into the phase location of 
BFCs in pulse profiles. The precursor location of both BFCs 
observed in J1012+5307 (at <j> ~ 71° and 245° in Fig. HJ 
is reasonable, because for fast-rotating pulsars with large 
dipole inclination, radiation emitted from any altitude can- 
not be detected far on the trailing side of the main pulse. Be- 
cause of the effects of aberration and retardation, the latest- 
possible phase of detection is equal to the caustic phase of 
16° s 2 / sin 2 a + 0/, where s is the footprint parameter of 
B-field lines and 4>f is the absolute fiducial phase (Dyks, 
Wright & Demorest 2009, in preparation). J0437-4715 does 
not obey this 'precursor rule', probably because of much 
smaller dipole inclination (a -C 90°, £ <C 90°). 



6.1.2 Azimuthal separation of orthogonal modes 

The plane of ip = 0, (orthogonal to the page in Fig. [7^.), 
is the plane of symmetry of the curvature beam. Such fan 
beam has therefore intrinsic mirror symmetry with respect 
to the plane of B-field lines. In the case of the partial contri- 
bution of the ||-mode (case ||B in Fig. [7}, the two orthogonal 
polarisation modes are naturally separated in magnetic az- 
imuth: the 1 1 -mode dominates in the plane of the B-field 
lines, whereas the _L-mode dominates on both sides of the 
parallel mode. This provides a very natural way to obtain 
the azimuthally-structured polarised carousels observed by 
Rankin et al. (2006) and Edwards et al. (2003). 

Let us assume the usual picture in which the polar re- 
gion is populated with thin plasma streams that are sepa- 
rated in magnetic azimuth and are anchored to the sparking 
spots (the carousel of Ruderman & Sutherland 1975). In 
such a case, the geometry of the curvature beam contami- 
nated by the ||-mode assumes the shape shown in Fig. 1141 
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(17) 



± mode 





Figure 14. In the case of the curvature radiation with partially- 
absorbed ||-mode, the carousel of electron streams (thick dashed) 
naturally produces polarisation pattern in which the two orthogo- 
nal modes dominate alternatingly in magnetic azimuth. The thick 
contours present the ||-mode, the thin contours are for the 
mode, and the B-field lines are dotted. The two horizontal lines 
mark two different ways in which our line of sight can leave the 
beam: in the case shown at the top the sightline cuts through the 
fixed-azimuth outer boundary of an azimuthally-limited emitter. 



The thin plasma streams (actually only their fragments that 
emit detectable radio waves) are marked with thick dashed 
lines. The ||-mode emission (dominating within, but not con- 
strained to, the planes of the streams) is delineated with 
thick solid contours. The _L-mode emission flanks the sides 
of the ||-mode beam (thin contours). The pattern naturally 
results from the intrinsic properties of the curvature beam 
with not-fully-absorbed |j-mode. Thus, the azimuthal sep- 
aration of modes can appear even with no refraction nor 
birefringence involved. 

Note that Fig. [14] presents the instantaneous geometry 
of emitting streams, and refers to single pulse data after they 
have been deconvolved into the image of polar carousel. In 
the averaged pulse profile, the structure of Fig. [14] is likely 
to become smeared out by plasma drift effects. The exact 
outcome of such averaging depends on the relative speed of 
drift motion and pulsar rotation. 



6.2 Edge depolarisation of profiles 

Rankin & Ramachandran (2003) emphasized the significant 
observation that the outer edges of pulsar profiles are al- 
most totally depolarised. Thus, some high-level conspiracy 
between the modes seems to be required. The curvature ra- 
diation is an extremely likely cause of this conspiracy, be- 
cause it intrinsically emits equal amounts of the orthogo- 
nal polarisation states in the outer edges of the curvature 
beam. For ip S> VWx, we have £ = I/7 2 + 4> 2 — ^ ■, an d 
Ki(y) — Ki(y), so the eq. @ implies: 



This ratio, as well as the polarisation degree L/I, are shown 
as a function of tp in Fig. [7p. The decrease of L/I with the 
angular distance from the symmetry plane (increasing 
is readily apparent. 

The well-known problem is that the emitted curvature 
radiation is polarised elliptically. It is not only limited to 
the two orthogonal modes. So some mechanism is needed to 
filter out the purely-orthogonal modes0 Whatever mecha- 
nism actually does it, a seemingly natural way to have equal 
amounts of two strictly-orthogonal modes, is to filter them 
proportionally out of a beam that already (intrinsically) in- 
cludes equal amounts of radiation in the two orthogonal po- 
larisation states. The curvature emission is then an ideal 
candidate, because in the outer wings of the elementary 
beam, it intrinsically produces equal amounts of radiation in 
the appropriate polarisation states. We therefore claim that 
the observed edge depolarisation has its origin in the intrin- 
sically equal emissivity in both polarisation states. Our line 
of sight can pass through the pulsar beam shown in Fig. [14] 
in two ways, marked with the horizontal lines. In the bottom 
case (grazing the edge of the carousel) it is probably a bit 
difficult to guess the averaged polarisation degree without 
numerical simulation. In the top case, however, the polari- 
sation degree at the edge of the profile will decrease, because 
our line of sight moves away from the fixed-azimuth plane 
of the last-seen emitter (|^| in Fig. [7p increases). 



6.3 The origin of double conal profiles 

An interesting question is what happens when the stream 
has noticeable azimuthal width, that is comparable to (or 
larger than) the opening angle ipm&x.- We argue that a cut 
through such stream produces the double conal profiles 
(type D in the classification scheme of Rankin 1983). A ge- 
ometry of such a cut is presented in the bottom right corner 
of Fig. 1181 The brightening of the edges of the profile is 
caused by the increased contribution of ||-mode. The near- 
absence of the 1 1 -mode in the inner parts of the profile leads 
to the smaller flux and large polarisation degree. The pres- 
ence of the ||-mode at the flanks of the profile also produces 
the depolarisation. In some cases, when sufficiently large 
amount of the ||-mode gets through, orthogonal jumps of 
averaged polarisation degree are observed in the outer wings 
of a pulse (eg. B1133+16, B2020+28, B2053+21, Hankins & 
Rankin 2008, hereafter HR08). 

The nearly-complete depolarisation of profile's edges is 
ensured by the fact that our line of sight passes through a re- 
gion with limited magnetic azimuth (as can be seen in Fig. 1181 
the sightline enters/exits the fixed-azimuth sides of the elon- 
gated wedge). It is not passing through a conal ring centered 
at the dipole axis. The brightest parts of the emission region 
(and the peaks of a profile) correspond to the fixed value of 
4>m = <j> l m on the leading edge and to <f> m = </>m on the trail- 
ing edge of the wedge-shaped stream. When our line of sight 



9 Several observations show that the separation can be very in- 
efficient, and the observed radiation is also polarised elliptically 
(eg. B2044+15, outer wings of B0525+21, B0751+32 and many 
other cases, see eg. Hankins & Rankin 2008; Edwards & Stappers 
2004). 
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departures from the fixed-ci m edge, it moves away from the 
center of the curvature beam shown in Fig. (7£i. The polari- 
sation degree is therefore decreasing as shown in Fig. [TJd. 

6.3.1 The origin of the radius-to-frequency mapping 

The interpretation here implies that the observed pulse 
width W , as measured between the profile peaks is fixed 
by the extent of the emitter/stream in the magnetic az- 
imuth: Ad> m — (j>m — (j) 1 ^- In the zeroth-order consideration 
this implies the lack of pulse broadening at low frequen- 
cies Some of the D-type pulsars indeed exhibit no widen- 
ing within a huge frequency range (eg. B0834+06 between 
~ 25 and 2370 MHz, see Fig. 4 in HR08), but some definitely 
do (eg. B0301+19). 

However, on a closer inspection the widening (mis- 
interpreted as RFM) can appear at low frequencies, since 
the plasma density in the stream falls off with the transverse 
distance from the stream center/axis. So we assume that the 
basic principle of the broadening is the same as usual: low- 
fobs radiation emerges from regions with lower density. How- 
ever, the lower density regions in our scenario are simply the 
more extreme outskirts of the stream, ie. the lower-density 
region located further away from the stream axis. Thus, we 
replace the 'radius-to-frequency' mapping with the 'stream- 
diameter-to-frequency' mapping. This interpretation implies 
that all radio frequencies come from the same (or similar) 
range of altitudes. 

6.3.2 The origin of the S-swing of polarisation angle 

The small width A<f> m of the emitter in the magnetic az- 
imuth implies that the projected direction of magnetic field 
does not change much within the pulse window. Thus, our 
scenario implies almost no change of polarisation angle (PA) 
due to the effect predicted by the rotating-vector-model 
(RVM) of Radhakrishnan & Cooke (1969). The PA curve 
is predicted to stay flat only if it is dominated by the RVM 
effect, ie. if the observed PA is determined by the projection 
of B on the sky. If the PA is governed by non-RVM factors, 
then behaviour that is in general different from the RVM 
curve is expected. 

The non-RVM behaviour is indeed observed in sev- 
eral pulsars with double peaks. For example, this is the 
case for B0834+26, B1919+21 (HR08), as well as for the 
geodetically-precessing pulsar B1913+16 (see fig. 7 in Weis- 
berg & Taylor 2002). However, the non-RVM origin seems 
to be plain inconsistent with the paramount examples of the 
RVM exhibited by two well-known representatives of the D- 
class: B0301+19 and B0525+21. They are widely believed 
to present textbook cases of the RVM-generated PA curve. 

We claim that here again we become victims of mali- 
ciously deceiving degeneration of Nature. The RVM model 
is really what the name says: a curve determined by the sky- 
projection of a vector that rotates under our line of sight. 
The projection of the circumpolar B-field lines is just one ex- 
ample of such a vector. The other examples, not even slightly 

10 In this paper we reserve the name "radius-to-frequency map- 
ping" (RFM) only to the geometric interpretation of the broad- 
ening. It is not used in reference to the observed phenomenon. 





A 


B 


C 


B 


A 


a) ~ 


- -5 


^ 1 ^ 



1 


5 b) 



-5 5 

pulse longitude <p [deg] 

Figure 15. a) Decomposition of a D-class profile into elementary 
beams emitted from different locations in the plasma stream. The 
_L-mode-dominated beam (thick line) is emitted from the center, 
and the ||-mode-contaminated, asymmetric beams from the out- 
side edges (thin lines). In the phase-interval C the variations of 
polarisation angle are determined by the intrinsic properties of 
the elementary beam (see text for more details), b) Plasma den- 
sity distribution sampled by our line of sight. The short bars mark 
the emission points for the curves in panel a. 



less natural than the projected B, are the vectors of electric 
field in the radiation pattern surrounding the instantaneous 
velocity of a relativistic particle. We argue that the famous 
S-swing of PA in the B0525+21 (and other D pulsars) is 
mostly caused by intrinsic effects. It is due to the cut of 
our sightline through the elementary (microscopic) _L-mode 
curvature beam. 

A rather complicated code would be needed to prove it 
numerically. Instead of that, we choose to delinate the mech- 
anism qualitatively. In Fig. 115b the observed pulse profile of 
a D-class pulsar is decomposed into the elementary beams, 
emitted from three locations in the stream. The thick curve 
presents the bifurcated _L-mode-dominated beam emitted 
from the center of the stream. The thin lines on the left- and 
right-hand side present the beams emitted from the outside 
left and right edge of the beam, respectively. The edge beams 
have their outside parts brightened by the presence of the ||- 
mode. Within the phase interval marked with the letter 'A', 
the polarisation properties are determined by relative pro- 
portions of the two modes. If the segregation of the modes 
is not efficient in 100%, a substantial elliptical polarisation 
appears there (as is often observed). In the phase interval 
'B' the flux observed at each phase always corresponds to 
the maximum of different elementary beams, emitted at dif- 
ferent locations in the stream. The changes of PA within the 
interval 'B' are therefore dominated by the changes of the B 
field direction. Thus, it is in the phase interval B, where the 
emission from different elementary beams is averaged out 
and can reveal the underlying geometry of magnetic field. 
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Figure 16. The origin of the PA swing in D-type radio pulsars. A 
head-on view of the _L-mode curvature beam with the orientation 
of polarisation vectors marked at selected places. The passage 
of observer's line of sight (horizontal line) leads to identical PA 
swing as in the case of the RVM model. 



In this interval the observed average PA curve can follow 
the RVM model. This can only happen if the width of the 
phase interval B is noticeably large (the width is determined 
by the relative width of the stream and the opening angle 
V'max of the _L-mode beam) . The corresponding parts of the 
PA curve of B0525+21 are indeed very flat (and highly po- 
larised because of the domination of the _L-mode) . 

In the phase interval 'C the situation changes com- 
pletely: for <f) < the flux at each phase is now dominated 
by radiation from a single elementary beam. Our line of sight 
does not sample different beams emitted from different loca- 
tions in the stream. The sightline moves through the inner 
parts of the elementary beam, and the changes of PA are 
determined by the changes of polarisation direction intrin- 
sic to the _L-mode beam (note that in this phase interval 
the received radiation is associated with a fixed direction 
of B). For (j> > (still in the C interval), the line of sight 
again samples a single elementary beam. Therefore, in the 
innermost parts of the double conal profiles, the PA is deter- 
mined by the intrinsic polarisation properties of the _L-mode 
curvature beam. As shown in Fig. 1161 the orientation of the 
polarisation vectors (vectors of the electric field E w of the 
emitted wave) in the instantaneous _L-mode-beam ensures 
the same variations of the PA as the projected magnetic field 
in the RVM model. This interpretation is based on the in- 
trinsic properties of the emitted radiation and assumes that 
they are not affected by any adiabatic walking below the 
polarisation-limiting radius (eg. Cheng & Ruderman 1979; 
Melrose 1979; Lyubarsky 2002). 

We thereby disclose a perverse situation, in which the 
central parts of the observed PA curve, that resemble the 
RVM shape so closely, have the origin that has nothing to 
do with the RVM idea. On the contrary, the outer flat parts 
of the observed PA curve (interval B) are more likely to 
reflect the underlying orientation of B and to comply with 
the RVM origin. Had this been the case, the difference of PA 
measured in the outer wings of the profile A (PA) would pro- 
vide information about the azimuthal width of the stream: 



A<f> m ~ A (PA) (provided the PA observed in the A and B 
region is not biased/dominated by intrinsic effects). 

The same conclusion can be reached in a different way, 
by considering the case of an infinitely-thin stream, and not- 
ing that the observed PA curve should change continously 
for small increase of the stream width. Specifically, for a 
stream of electrons flowing only along a single B-field line, 
the observed PA curve definitely has to be determined fully 
by the intrinsic swing shown in Fig. 1161 This is because the 
spatial extent of the emitter is a delta function and there 
is no variation of the underlying direction of B. A small 
increase of the diameter of the stream must, by continu- 
ity, produce only small changes of the initial PA curve. The 
increasing spread of B within the stream of increasing thick- 
ness will therefore produce only a gradual change in the ob- 
served PA curve. We therefore conclude that in general the 
observed PA curves of pulsars are a convolution of intrinsic 
PA and B-field projection. 

Since the opening angle of the _L-modecurvature beam 
decreases with increasing frequency, the PA curve is ex- 
pected to depend on u obB , which is indeed the case 
for B0525+21 (and other objects). Because the _L-mode- 
radiation has a minimum in the profile center (Fig. I15p . any 
remnants of the ||-mode emission can produce a minimum of 
polarisation degree there. Even such detail is also observed 
in B0525+21 (eg. fig. 1 in HR08). In a more general case 
(than the primitive picture of Fig. I15[) . a large variety of 
PA curves is expected, depending on the 3D distribution of 
plasma density in the stream, on the geometry of passage 
of our sightline through the beam, on the ratio of both or- 
thogonal modes, and on the size of the elementary curvature 
beam. For realistic density distribution the phase intervals 
B and C can overlap in phase, in which case the PA recorded 
at a given phase can be partly determined by the intrinsic 
PA and partly by the projection of B. The subject deserves 
numerical study. 

Interesting examples of intrinsic variations of PA curve 
have also been observed under core components, where the 
polarisation has been found to depend on the subpulse 
brightness (B1237+25, Srostlik & Rankin 2005; B0329+54, 
Mitra, Rankin & Gupta 2007). The core emission has been 
found to be polarised orthogonally with respect to the outer 
components (Srostlik & Rankin 2005) which suggests ||- 
mode origin of the core (Gil & Snakowski 1990). 

We conclude that the observed PA curve (at least for 
D-type pulsars) is strongly influenced by the intrinsic po- 
larisation properties of the _L-mode-beam. Even in the most 
RVM-like looking fragments it does not have the RVM orgin. 
Though this interpretation invalidates the determination of 
the general geometric parameters (a, £) from the PA curves, 
it opens a way to interpret their ubiquitous distortions from 
the RVM-shape. 

6.3.3 Profile width and the geodetic precession 

Precessing pulsars give us the unique opportunity to view 
their emission beams in two dimensions. Especially pre- 
cious are the pulsars in relativistic binary systems, such 
as J0737-3039A/B, B1913+16, J1141-6545, B1534+12, or 
J1906+0746. They should exhibit the geodetic precession, 
the rate of which is well determined by the general rela- 
tivity. For the traditional, axially-symmetric pulsar beams, 
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Figure 17. Variations of pulse width W as a function of the 
viewing angle f in the stream model (solid lines; numbers give 
the value of <j>m = A0 m /2). The outermost two curves present 
the thinnest streams of width A0 m = 5° and 10° in the magnetic 
azimuth. The streams are centered at the main meridian and 
a = 60°. Dashed line presents the standard case of a conal beam 
with the opening angle 28 m = 10° . 

such precession should generally produce clear changes of 
profile width. Several detailed studies of the phenomenon 
(Kramer 1998; Clifton & Weisberg 2008; Lorimer et al. 2006; 
Burgay et al. 2005) have reported that the profiles refuse to 
change at the expected rate. This suspicious observation can 
be explained in several ways, eg. by near-alignment of the 
spin axis and orbital angular momentum, or by the acci- 
dental approximate co-planarity of the rotation axis, orbital 
momentum, and line of sight (see Clifton & Weisberg 2008). 
However, none of the explanations seems natural or likely, 
especially in view of the increasing ubiquity of the problem. 

We argue that the real reason for the profile stability is 
the completely different topology of pulsar beams: our sight- 
line is cutting through the split/structured fan beams gen- 
erated by elongated plasma streams, not through the narrow 
cones. The pulse profile of 131913+16 indeed has the D-class 
shape and exhibits strongly distorted (non-RVM) PA curve 
(fig. 7 in Weisberg & Taylor 2002). However, our fan-beam 
interpretation is likely to hold not only for the D-class pro- 
files. The pulsar J1906+0746 (Lorimer et al. 2006) belongs 
to a single-component class (S type), and exhibits interest- 
ing profile evolution: despite its interpulse clearly appearing 
within a few-years-long observation time, the width of the 
main pulse does not vary considerably. A cut through fan- 
beams explains such behaviour naturally. 

In the traditional conal model the observed pulse width 
W = 20 is modelled by assuming the fixed opening angle 
8 m of the conal beam in the equation: 

, cos 9 m - cos a cos ( . 

cos0= : : — . (18) 

sm a sin (" 

In the fan beam scenario it is the azimuth 4> m of a wedge- 



shaped stream that needs to be fixed, whereas 8 m in the 
above formula changes with time, complying with the equa- 
tion: 

cos £ = cos(7r — <j>m) sin a sin 9 m + cos a cos 8 m . (19) 

The assumption of W — 2(j> is only correct for a stream that 
is symmetric with respect to the main meridian (containing 
the rotation and magnetic axes). In general, <^ 7^ and 
the width is equal to W = <f>(4>m) ~ 0(0™,) i with <f> given 
by eq. (|18p . The instantaneous value of cos# m needs to be 
calculated from the quadratic equation (|19[1 and used in (|18l) . 
This leads to the much slower variations of W with (, as 
can be seen in Fig. [T7] in which the solid lines show W for 
streams of various azimuthal width 8cj> m = 2(j> m and a — 
60°. The quasi-linear changes of W , characteristic for the 
observations of geodetically-precessing pulsars, are readily 
apparent for narrow streams (outermost lines; see the values 
of (j)m marked at each line). The changes of W for a conal 
beam with 8 m — 10° are shown for comparison with the 
dashed line. 

Therefore, because of the geodetic precession we 
have already traversed a long way through the beam of 
B1913+16. The profile has not evolved much because we 
cut through a fan-shaped beam. We argue that the sta- 
bility of profile in both components of the double pulsar 
J0737— 3039A/B has the same origin, despite the profile clas- 
sification for the A star is far from clear. 



6.4 Sketch of a new classification scheme 

6.4-1 Neither conal nor patchy: spoke-like, wedge-like, and 
stream-like 

We have come to the conclusion that bifurcated compo- 
nents in radio-pulsar profiles are produced by sightline cuts 
through fan-beams of radiation from elongated streams. 
There need to be stream-like emitters in the pulsar magne- 
tosphere, as shown in the bottom-left corner of Fig. 1181 An 
important question is what is the length of the streams and 
what is the range of magnetic colatitudes that the fan beams 
subtend. At least in some cases the fan-beams must extend 
quite far from the dipole axis, to explain the large separation 
of bifurcated components from the main parts of the profile 
(main pulse, interpulse). The large extent could also explain 
the abundance of pulsars with interpulses (eg. Weltevrede 
& Johnston 2008b). If the streams are not exactly centered 
at the main meridian, the main-pulse-interpulse separation 
can differ from 180° as it sometimes happens to be observed 
(eg. in B1055-52, Weltevrede & Wright 2009). 

The components of the BFC class are observed at 
frequency-independent location in pulse profiles and widen 
with increasing u h B roughly according to Abf c oc . 

For the price of rejecting the RVM wisdom, we have 
also come to the conclusion, that the double conal profiles 
(class D) are also created by a cut through azimuthally- 
limited emitter, that is a bit broader than in the BFC case 
(right bottom corner of Fig. I18|) . The width of such pro- 
files is mostly determined by the fixed width of the stream. 
Any low-frequency widening (of non-RFM nature) can re- 
sult from density gradient in the outskirts of the stream. The 
intrinsic u~^ 3 widening can also contribute to the apparent 
effect. 
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Figure 18. A head-on view down the 
magnetic pole, presenting the geome- 
try of possible emission regions. Some 
of them are plotted at two frequen- 
cies: a lower one (dashed contours) 
and a higher one (solid contours). The 
standard conal view of profile gener- 
ation (middle sightline) is contrasted 
with the stream origin that we ad- 
vance in this paper (bottom sight- 
lines). The stream-cut origin is likely 
not restricted to the D-class pulsars, 
as suggested in the upper part of the 
figure. 



The fan-beam origin of pulse profiles is probably not 
limited to the D class. The profile of B1929+10 is likely cre- 
ated by a cut through a structured stream, or through a few 
nearby streams. At 430 MHz, its profile has a pronounced 
emission component on the trailing side of the maximum, 
which separates from the maximum at a rate roughly con- 
sistent with f obs . Moreover, in the minimum between the 
component and the central maximum, noticeable amount 
of orthogonal polarisation mode appears (Rankin & Rath- 
nasree 1997). This kind of corellation (appearance of an- 
other polarisation mode at a minimum in total flux) could 
in principle be produced by any process that removes the 
dominating mode. However, here we find again that the 
_L-mode-dominated curvature beam (case ||B in Fig. [7^) 
has the appropriate intrinsic property to produce the phe- 
nomenon (domination of the ||-mode at the minimum of 
the ±-mode)EI We therefore suggest that the 'multiple' ap- 
pearance of the profile of B 1929+ 10 is partially caused by 
the intrinsic doubleness of the curvature beam, and the real 
number of underlying streams is smaller than the number of 
observed components. If the radiation in the pedestal emis- 
sion is averaged out to reveal the RVM polarisation angle, 
then this pedestal PA curve should be discordant with the 
PA curve observed under the main pulse (because the lat- 
ter is dominated by intrinsic effects). This likely explains 
the difficulty in fitting the full-period interval with a single 
PA curve (Rankin & Rathnasree 1997; Everett & Weisberg 
2001). B0950+08 may belong to the same, B1929+10-like 
class: its profile tends to bifurcate at low v bs- The extended 
leading-bridge of this profile is likely caused by a stream that 
extends roughly tangentially to the trajectory made by our 
line of sight. The bridge culminates in an interpulse which 
again reveals doubleness (McLaughlin & Rankin 2004). 

A paramount example of stream-cut-generated profile 



11 Another beautiful example of such effect is seen in B1745— 12 
(fig. 7 in Mitra, Sarala, & Rankin 2004). 



is provided by J0437— 4715 in which the tendency to 'look 
double' is revealed both in the complicated central parts and 
in the conal components. Whereas the core can be produced 
in the usual way (a cut through the dipole-axis-centered 
rings) it may also be created by streams flowing near the 
main meridian. The outer "conal" components are definitely 
produced by streams with large 4>m- Deconvolution of this 
(and other) profiles into the underlying emitting structures 
(streams) can only be done if the bifurcated nature of the 
elementary curvature beam is taken into account. 



7 SUMMARY 

The bifurcated component of J1012+5307 can be 
comparably-well fitted with both the parallel-acceleration 
emission and the curvature radiation. However, only the 
curvature radiation can match the data simultaneously at 
two frequencies. The two-planar topology of the _L-mode- 
dominated curvature beam makes it very easy to reproduce 
the observed large depth of double notches. It also explains 
the origin of bifurcated components in terms of very natural 
geometry (sightline cuts through plasma streams). This 
leads us to the firm conclusion that the observed coherent 
radio emission from pulsars is the curvature radiation. 

Components in radio pulsar profiles become wider, and 
finally bifurcate at low-frequencies because of the intrinsic 
behaviour of the curvature beam. The scale of the bifurca- 
tion makes it possible to determine the radius of curvature 
of electron trajectory in the emission region, provided global 
geometry parameters are known. 

Double "conal" profiles (D class) are not conal in the 
usual sense. They are produced by the sightline cut through 
slightly broader streams. 

The PA curves of such profiles are not determined by 
the projection of the _B-field on the sky plane. Even when 
they resemble the RVM shape (which they do rarely), the PA 
curves are mostly determined by the microphysical polarisa- 
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tion structure of the _L-mode curvature pattern. The S-swing 
created by a cut through such single pattern is exactly the 
same as in the RVM case. Angular/spatial convolution of 
more patterns leads to a variety of observed PA curves. 

The low-frequency broadening of D-type profiles results 
from plasma density gradient at the outer surface of the 
stream, and from the intrinsic broadening of the curvature 
beam. It does not reflect the emission radius (it does not 
have the RFM origin). 

The narrow-stream geometry of the emitter implies fan- 
like shapes of pulsar emission beams. This naturally explains 
the extremely slow changes of pulse width due to the geode- 
tic precession. 

The existence of fan-shaped beams increases the like- 
lihood of observing an interpulse. It also implies that the 
number of pulsars in the Galaxy may be smaller than esti- 
mated with the conal beams. 

The results of this paper invalidate the traditional 
methods of determination of the radio emission altitude at 
least for D-class pulsars. This refers both to the geometrical 
method, and the one based on the relativistic-shift of PA 
curve. 

The curvature-radiation beam with the partially- 
attenuated ||-mode has several features that are qualita- 
tively consistent with many observed properties of pulsars. 
The features include the non-zero width of the curvature 
beam, the mirror symmetry with respect to the B-field line 
plane, the 'striped' (aziumuthally-separated) polarisation 
structure, and the intrinsic 1:1 ratio of the orthogonal polar- 
isation modes at the beam edge. They enable us to roughly 
comprehend the multiple imaging, azimuthal separation of 
modes, edge depolarisation, and the polarisation structure 
of D-type profiles. 

The stream origin is likely for other types of pulsar pro- 
files. In attempts to decipher specific cases, it may be more 
efficient to think in terms of spoke-like, or wedge-like emit- 
ter, not in terms of 'conal or patchy'. 



8 DISCUSSION 

Several coherency mechanisms for the amplification of the 
curvature radiation had been proposed. They are usually di- 
vided into the bunching processes, and the maser processes, 
and both are reportedly problematic (Melrose 2006). 

The elementary beam of the coherent curvature radi- 
ation, that we have invoked from the fits of Section [3] and 
from the pulsar properties of Section [5] suggests that the 
coherency is not due to the curvature-drift-driven maser 
(Zheleznyakov & Shaposhnikov 1979; Luo & Melrose 1992, 
hereafter LM92) nor due to the maser induced by the field- 
line torsion (Luo & Melrose 1995). It is true that the first 
model predicts the required two-planar topology of the beam 
and the correct frequency dependence of Abf c - However, the 
maser implies intrinsically asymmetric beam which does not 
seem to be consistent with the wonderful symmetry of the 
double features. This can be seen in fig. 1 of LM92. The am- 
plified parts of the non-coherent beam (determined by the 
negative absorption coefficient P in the bottom-right panel) 
are located asymmetrically with respect to the plane of elec- 
tron trajectory (at the angles marked 9\ and 82 in the fig- 
ure of LM92). The maser operating thanks to the B-field 



line torsion (Luo & Melrose 1995) amplifies the radiation 
at extremely small angles with respect to B. The bunching- 
induced coherency (eg. due to the two-stream instability, 
Ruderman & Sutherland 1975) seems to have less problems 
with the quasi-isotropic amplification of the non-coherent 
beam. Recently this idea was studied by Gil, Lyubarsky & 
Melikidze (2004) in their version of 'soliton bunching'. In 
their figure 3 they present elementary beams of the ±-mode 
radiation that are fairly similar to the BFC of J1012+5307. 
It is less clear if their mechanism can amplify comparable 
amounts of the ||-mode in the outer parts of the beam. 

The idea of "fan" beams appeared from time to time 
in pulsar literature. Eg. Narayan & Vivekanand (1983) con- 
sidered an elongated beam to model the notorious profile 
of B0950+08. However, the beam elongation in their model 
(Fig. 2 therein) was still caused by a deformation (squeezing) 
of a conal beam and did not have the stream-like topology. 
It is peculiar that in the past the fan-beam origin of profile 
components has not been emphasized sufficiently enough to 
spark development of a spoke-like model of pulsar beams. 
The intrinsic fan-beam topology of the curvature radiation 
has been known from the early days of synchrotron radiation 
theory, and the tendency of curvature radiation to produce 
double components has even been incidentally illustrated in 
the literature (eg. figs. 2 and 3 in Ahmadi & Gangadhara 
2002). 

We conclude that the long-sought Rosetta Stone needed 
to decipher the nature of pulsar radio emission has finally 
been identified as the double features in averaged pulse pro- 
files. 
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